Magnetic order in the quasi-two-dimensional easy-plane XXZ model 
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A Green's-function theory of antiferromagnetic short-range and long-range order (LRO) in the 
S — 1/2 quasi- two-dimensional easy-plane XXZ model is presented. As the main new result, two 
phase transitions due to the combined influence of spatial and spin anisotropy are found, where below 
the higher and lower Neel temperature there occurs LRO in the transverse and in both the transverse 
and longitudinal spin correlators, respectively. Comparing the theory with neutron-scattering data 
for the correlation length of La2Cu04, a very good agreement in the whole temperature dependence 
is obtained. Moreover, for La2Cu04, S^CuChCh, and Cao.8sSro.i5Cu02 the second phase with 
longitudinal LRO is predicted to appear far below room temperature. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.-s 



I. INTRODUCTION 

To understand the unconventional behavior of high-T c 
superconductors, which is mainly ascribed to a strong 
antiferromagnetic (AFM) short-range order (SRO), 
the magnetic properties of the quasi-two-dimensional 
(2D) parent compounds were probed preferably by 
neutron scattering and NMR experiments, e.g., on 
La 2 Cu0 4 [1-4], Ca . 85 Sr . 15 CuO 2 [5], YBa 2 Cu 3 6+x 
(x < 0.4) [6] and L 2 Cu0 4 (L=Nd, Pr) [7]. In partic- 
ular, the staggered magnetization [1-3,5,6] and the AFM 
correlation length [2,4] were investigated. To analyze 
the experimental data, quasi-2D Heisenberg models in- 
cluding a weak spin anisotropy (easy-plane XXZ models) 
were considered [1-7] and treated by linear spin- wave the- 
ory [3,5,6] and a Schwinger-boson approach [1]. Recently, 
a quasi-2D anisotropic (easy-axis) Heisenberg model was 
studied within a boson- fermion mean-field theory [8]. 
However, those auxiliary-field approaches and most of 
the spin-wave theories for related spin models (for refer- 
ences to most of the early work, see Ref. [9]) yield reason- 
able results only at sufficiently low temperatures, since 
the temperature-dependent SRO is not adequately taken 
into account. To provide a good description of magnetic 
SRO at arbitrary temperatures, a Green's-function the- 
ory based on the projection method was developed for 
the spatially anisotropic Heisenberg models [9-11] and 
for the ID [12] and 2D easy-plane XXZ models [13] . 

In this paper we extend our previous work to the S = 
1/2 quasi-2D easy-plane XXZ model 
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and (i,j) z denote nearest-neighbor (NN) 
sites in the xy plane and along the z direction of a simple 
cubic lattice, respectively, R z — J z /J < 1 describes the 
AFM interplane coupling (throughout we set J = 1), 
and the spin anisotropy parameter A is considered in the 
AFM region < A < 1. In the limiting cases A = 1 and 
R z = the model (1) reduces to the models studied in 
Refs. [9] and [13], respectively. 

Our Green's-function theory outlined in the Appendix 
is based on an approximate time evolution of spin oper- 
ators, e.g., — 5q = (Wq Z ) 2 5q, resulting from a decoupling 
of three-spin operator products which is improved by the 
introduction of vertex parameters. For the temperature 
dependence of some vertex parameters assumptions are 
made which are detailed and motivated in the Appendix. 
Within this theory we examine the combined effects of 
spatial and spin anisotropy on the AFM long-range order 
(LRO) at T = (Sec. II), the Neel transition tempera- 
ture, and on the AFM correlation length (Sec. III). In 
Sec. IV we compare our results with experiments. For 
La 2 Cu04, Sr 2 Cu0 2 Cl 2 , and Ca .85Sr .i5CuO 2 our the- 
ory predicts, in addition to the usual Neel transition, 
a further transition far below room temperature, where 
the spin correlators between the z-components start to 
develop AFM LRO. 



II. GROUND-STATE LONG-RANGE ORDER 

Analyzing the magnetic LRO described by the stag- 
gered magnetizations m u [v = ±,zz; cf. Eq. (11)] at 
T = as functions of R z and A, we obtain transverse 
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LRO in the whole parameter region considered and two 
solutions differing in the existence of longitudinal LRO. 
That is, we obtain a phase with m zz — (phase I) and a 
phase with m zz ^ (phase II), where in both phases we 
have m H 7^ 0. The stabilization of longitudinal LRO 
in the easy-plane region by the interplane coupling may 
be due to the reduction of quantum spin fluctuations in 
higher dimensions. 

In Fig. 1 the R z — A phase diagram is shown, where for 
R z 7^ the transition across the phase boundary denoted 
respectively by i? ZjC (A) and A C (R Z ) is found to be of sec- 
ond order (cf. inset). Let us consider the phase transi- 
tion in the vicinity of the critical point (A,R Z ) = (1,0) 
in more detail. In our approach the solution for m zz 
turns out to depend sensitively on the input data for 
<9e(A, \)/dA = e' used to determine the vertex param- 
eter af z , where e(A, 1) denotes the ground-state energy 
of the 2D XXZ model (see Appendix). Taking for e(A, 1) 
the exact diagonalization (ED) data on lattices with up 
to 36 sites (without finite-size scaling) from Rcf. [13], we 
obtain lim^^o A C (R Z ) = A — 0.958. On the other 
hand, taking the Monte Carlo (MC) data from Ref. [14], 
which have to be interpolated between the few available 
points A = 0, ±0.5, ±1, we get limA^i i? z , c (A) = i? Zi o = 
4.08 x 10~ 2 . However, as is well known, at R z = there 
is no longitudinal LRO for < A < 1 which means that 
we must have Aq = 1. Moreover, we expect the inter- 
plane coupling to be a relevant perturbation with respect 
to the stabilization of LRO analogous to the situation 
in the 2D spatially anisotropic Heisenberg model [11,15], 
where LRO at a finite arbitrary small interchain coupling 
was found [15]. Therefore, we make the reasonable as- 
sumption R z , = 0. To fulfill the requirements A = 1 
and R z .o — simultaneously, so that the phase boundary 
touches the critical point (A,i? z ) = (1,0), as input for 
e' we use a linear combination of the ED [13] and Monte 
Carlo data [14], e' = xe' ED + (1 — x)e' MC . We find that 
the above requirements may be fulfilled, if x is chosen as 
x=0.44. 

In the limit A = 1 we have rotational symmetry 
(C zz = C+-/2), so that V2m zz (R z ) = m +L (R z ) = 
y // 2/3m(R z ) with m defined as in Refs. [9-11]. In this 
limit, our result for m H (R z ) agrees with that of Ref. [9]. 
As can be seen in the inset, the effects of spin anisotropy 
on the longitudinal and transverse LRO are opposite: 
We have dm zz /dA > 0, whereas <9to h /dA < which 
agrees, at R z — 0, with the Monte Carlo data [14] and 
the results of Rcf. [13]. 



III. FINITE-TEMPERATURE PROPERTIES 
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FIG. 1.: R z — A phase diagram and transverse and lon- 
gitudinal zero-temperature magnetizations (inset) in the 
quasi-2D easy-plane XXZ model. 



order phase transitions at (R Z ,A) and T^ Z (R Z , A) 
[m v {T v N ) = 0] with > T Z N Z and T Z N Z {R Z . C ) = 0. 

Figure 2 shows the Neel temperatures as functions of 
R Zl i.e., the T — R z phase diagram for different spin 
anisotropics. For R z — we obtain T£~ = 0, in agree- 
ment with the Mermin- Wagner theorem. At A = 1 we 
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Considering the AFM LRO in the phases I and II 
\m v (T = 0) 7^ at R z > R z , c ] at nonzero temperatures, 
the solution of the self-consistency equations (5), sup- 
plemented by the conditions for the vertex parameters 
[cf. Eqs. (12) and (15) to (17)], results in two second- 



FIG. 2.: T — R z phase diagram in the quasi-2D easy- 
plane XXZ model. Below the Neel temperature T^ z the 
phase with longitudinal long-range order becomes stable. 
The curves below the solid line belong to T^ z . 
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Tn. As compared with the results 
of Ref. [9], our values for TV are somewhat higher (by 
about 9%) due to numerical uncertainties. On the other 
hand, in comparison with previous RPA and mean-field 
approaches (cf. Ref. [9]) our Neel temperatures are re- 
duced by the improved description of SRO. For example, 
the Tjv values found by the Schwinger-boson approach of 
Ref. [16] exceed our results by a factor of about 1.7 on 
the average. 

Concerning the influence of spin anisotropy on the 
Neel transitions, we obtain dT^~/dA < and 
> 0, corresponding to the A dependence 
of m v (cf. inset of Fig. 1). The dependence on 
A of is in qualitative agreement with the be- 

havior found in previous approaches [1,2]. There, 
Tjv (being identified with T£~) is given as T/v/J = 
-27rM {ln|4a e// /[7r 2 M ln(4a e// /7r)]|}- 1 (Ref. [1]) with 
Mo = 0.3, a e ff — Aa xy + 2R Zl and a xy = 1 — A or as 
Tn/J = — 47T/9s(lna e //)~ 1 (Ref. [2]), where p s is the 
spin stiffness. However, in contrast to those mean-field 
(Schwinger boson) results, in our theory the combined in- 
fluence of spatial and spin anisotropy on cannot be 
expressed in terms of a single effective parameter. Con- 
sidering the variations SR Z and 5a xy under the condition 
ST^~ = we get sgn(5a xy ) = — sgn(SR z ). Whereas the 
Tjv formulas quoted above yield Sa xy — —5R z /2 for all 
R Zl from Fig. 2 we obtain a R z dependent relation be- 
tween Sa xy and SR Z . 

In Fig. 3, our numerical results for the temperature 
dependence of the magnetization vn v are depicted. They 
may be described by a T 2 decrease at low enough temper- 



atures (3D behavior) and by m u (T) oc (1 - T/T^) 1 / 2 for 
temperatures close to T^, as was also found in the A = 1 
limit [9] . The influence of spatial and spin anisotropy on 
m u {T) is analogous to that on m l/ (T — 0) shown in the 
inset of Fig. 1. 

In the inset of Fig. 3 the inverse AFM intraplane cor- 
relation lengths above are plotted, where the effects 
of the intcrplanc coupling and spin anisotropy are visi- 
ble. In the vicinity of the temperature dependence 
of (CnfjT) -1 changes from an exponential law in the 2D 
case (T^~ = 0) to a linear behavior for R z > 0. Equally, 
(£ccj/) _1 near T^f behaves as T—T^ z . According to the in- 
fluence of spin anisotropy on Tjy we get d^ xy / dA > and 
d£+ y ~ /dA < 0. The behavior of qualitatively agrees 
with the anisotropy dependence of the mean-field expres- 
sion for £ (being identified with £, x ~) given in Ref. [2], 

£ = £o(l — a e//Co) -1 ^ 2 ' where £o denotes the correlation 
length for R z — and A = 1. 

Finally let us consider the heuristic relation be- 
tween T£~ and the transverse 2D correlation length 
at Tjy - which is often used in describing the exper- 
imental data [17] and is given by Q(R Z ,A = 1) = 
0.25 with Q(R Z ,A) = R z (T+-)- 1 [m+-(T = 0, R z = 
0,A)£+-(T = T+~,R Z = 0,A)] 2 . By our results, 
Q(R Z ,A) in the experimentally relevant region (cf. 
Sec. IV) 2 x 10~ 4 < R z <2x 10~ 2 and at A = 1(0.8) 
is found to vary between 0.21 (0.31) and 0.095 (0.23). 
That is, the heuristic estimate is roughly confirmed by 
our theory. 
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FIG. 3.: Staggered magnetizations versus temperature. 
The inset shows the inverse intraplane correlation lengths 
above the corresponding Neel temperatures. 



IV. COMPARISON WITH EXPERIMENTS 

Let us first compare our results for the transverse 
intraplane correlation length with the neutron- 

scattering data on La 2 Cu04 (Ref. [4]) in the range 
340 K< T <820 K plotted in Fig. 4. Based on the 2D 
Heisenberg model (A = 1), in Ref. [18] the exchange en- 
ergy J was determined by a least-squares fit (a=3.79 A), 
where for our choice of the vertex parameters the re- 
alistic value J = 117 meV was found. Here, we fix 
this value and consider the effects of spatial and spin 
anisotropy on £ xy ~(T). The deviation of the theory for 
R z = and T < 550 K from the experimental data may 
be reduced by the inclusion of the interplane coupling, 
since £ xy (T N ) = 0. For A = 1 and T N = 325 K [1,4] 
we obtain R z = 1.2 x 10~ 3 , and the theoretical low- 
temperature £~y curve lies only somewhat above the ex- 
periments (cf. Fig. 4). Taking into account the spin 
anisotropy a xy — 1.5 x 10~ 4 [2] or a xy = 5.7 x 10~ 4 [4], 
for T+~ = 325 K we get R z = 3.0 x 10~ 4 . For those pa- 
rameters we obtain an excellent agreement between the- 
ory and experiment over the whole temperature region. 
Note that the theoretical curves for a xy = 2x 10~ 4 (cf. 
Fig. 4) and a xy = 2 x 10~ 3 agree within the accuracy of 
drawing. 
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Concerning our prediction of phase II with longi- 
tudinal LRO in La2CuC>4, for a xy = 1.5. x 1CP 4 
(5.7 x 1CP 4 ) we obtain the longitudinal zero-temperature 
magnetic moment [i zz = 2^m zz = 6.6 x 10~ 2 ^s 

(6.1 x 1CP 2 (U_b) as compared with the transverse moment 

+- 



we find T% z = 

rpzz 

N 



= 0.55^b- For the Neel temperature 
2.6 x 1(T 2 J (2.4 x lO- 2 J) equally to 
= 35 K (33 K) [J = 117 meV]. With regard to the 
experimental verification of the two phases, the magni- 
tude of the longitudinal moment (yu zz > 0.1^i H ) may be 
large enough to allow a separation between [i zz and ^ 
by polarized neutron-scattering studies on single crystals 
of La 2 Cu0 4 [19]. 

Next we consider the compound S^CuC^C^ which 
is the best experimental realization of an S = 1/2 2D 
Heisenberg antifcrromagnct, where J = 125 ± 6 meV, 
a xy = 1.4xl0- 4 , andT+~ = 256.5 K [20]. Keeping T+~ 
and a xy fixed (we take a xy = 2x 10~ 4 , as for La2Cu0 4 ) 
and choosing J = 125 meV, 120 meV, and 110 meV, 
we obtain the interplanc coupling R z = 2.0 x 10~ 5 , 
4.0 x 10~ 5 , and 1.0 x 10~ 4 , respectively. On the aver- 
age we have R z ~ 5 x 10~ 5 <C a xy in qualitative agree- 
ment with the estimate given in Ref. [20]. In Fig. 5 our 
results for the transverse intraplane correlation length 
(a = 3.967 A), where (C^)" 1 oc T-T+" near the transi- 
tion to phase I, are compared with the neutron-scattering 
data [20]. For J = 125 meV we obtain a very good agree- 
ment between theory and experiment at low enough tem- 
peratures [T < 400 K), whereas for J = 110 meV the 
agreement is good at higher temperatures. 

The data on the predicted phase II in Sr2Cu02Cl2 cal- 
culated for J = 125 meV, 120 meV, and 110 meV are ob- 
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tained as fi zz /^ B = 2.8 x 10~ 2 , 3.9 x 10~ 2 , and 5.5 x 10~ 2 
(for comparison, /x" 1 j hb = 0.54) and as T^ z = 7 K, 
12 K, and 18 K, respectively. As in La2Cu04, the mag- 
nitude of the longitudinal moment (/i zz ~ 4 x 10~ 2 /is ~ 
0.08^i H ) may be large enough to be detected by polar- 
ized neutron-scattering experiments on Sr2Cu02Cl2 sin- 
gle crystals. 
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FIG. 5.: Inverse antifcrromagnctic transverse intra- 
plane correlation length in S^CuO^C^ obtained by the 
neutron-scattering experiments of Ref. [20] (symbols) 
and from the theory for different exchange energies J. 
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Considering Cao.ssSro.isCuO^ (TV 
125 meV [5]), for A = 1 we obtain R z = 9.0 x 10"" as 
compared with R z ~ 2.5 x 10~ 2 resulting from a fit of 
the magnetization data [5] . Taking OL X y cls for La2Cu04, 
a xy = 1.5 x 10~ 4 [5], we get R z = 5.0 x 10~ 3 . For the 
zero-temperature magnetic moments and the Neel tem- 
perature we obtain [i zz = 0.16^b, M H = 0.57yUs and 
T^f — 190 K, respectively. Note that the longitudinal 
moment is more than twice as large (^ zz = 0.28//+ - ) 
as in La 2 Cu0 4 . However, contrary to La 2 Cu0 4 and 
Sr 2 Cu02Cl2, for Ca(Sr)Cu02 single crystals are not 
available, so that neutron-scattering data do not exist 
until now [19]. 
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FIG. 4.: Inverse antiferromagnetic transverse intraplane 
correlation length in La2Cu04 obtained by the neutron- 
scattering experiments of Ref. [4] (symbols) and from the 
theory for different spatial and spin anisotropies. 
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V. SUMMARY 

In this paper we presented a Green's-function theory 
for the quasi-2D easy-plane XXZ model allowing the cal- 
culation of all static magnetic properties at arbitrary 
temperatures, where we focused on the effects of spatial 
and spin anisotropy on the AFM LRO and the correla- 
tion length. As a qualitatively new result, for appropri- 
ate model parameters we obtained two phase transitions, 
where the paramagnetic phase with pronounced AFM 
SRO becomes unstable against a phase with transverse 
LRO only and, at a lower temperature, a phase with both 
transverse and longitudinal LRO. Comparing the the- 
ory with neutron-scattering experiments on the correla- 
tion length of La2Cu04, an excellent agreement is found. 
Furthermore, the second Neel transition (to the phase 
with longitudinal LRO) in La2Cu04, Sr 2 Cu0 2 Cl 2 , and 
Cao.85Sro.i5Cu02 is predicted to occur at about 30 K, 
10 K, and 190 K, respectively. Our goal is to stimulate a 
wider discussion and new experiments in this direction. 

Acknowledgments. The authors are greatly indebted 
to B. Kcimcr and R. Hayn for stimulating discussions. 



APPENDIX: THEORY OF SPIN 
SUSCEPTIBILITY 

The spin susceptibilities x+~(o>) = — ST q )) w an d 
X q » = Si q )) U) (((...;.. .» w denotes the two- 

time retarded commutator Green's function) are deter- 
mined by the projection method taking, as for the XXZ 
chain [12], the basis iS+) and (5 q , iS^), respectively. 
Wc obtain 



XqM = — 2 — ruyi ; 

H UJ Z — (Wq) 



v = H — , zz, 



(2) 



with the first spectral moments M+ = ([iS£ , SLq]) and 
Mq Z = {[iSq, Stq\) given by the exact expressions 

M+- = -4[C+ - (l - A 7q ) + 2C^ 0i0 (A - 7q )] (3) 

-2i? z [C + - 1 (l - Acos 9z ) + 2Qf M (A - cosq z )] 
M? = -4(7+ - (l - 7q ) - 2R Z C+^(1 - co S q z ) , (4) 



= (SqS t ), C zz — (SqS z ), r — 
ne x + me y + le z , and 7q = (cos q x + cos q y )/2. The spin 
correlators are calculated as 



1 M" r 



(5) 



where p(w q ) = (e< /T - The NN correlation func- 

tions are related to the internal energy per site e = 

2(Ci,o,o + ACfo,o) + -^z(Co,o,i + ACgQ ;1 ). 

The spectra o> q are calculated in the approximations 

-S+ = {oj+-) 2 S+ and -S* = (w$*) 2 5*, where prod- 
ucts of three spin operators in —Sf and — Sf along 



NN sequences are decoupled. Introducing ver- 

tex parameters in the spirit of the scheme by Shimahara 
and Takada [21] and extending the decouplings given in 
Rcfs. [11,9,12], we have 



S+S+S- = a+- lz (S+S-)S+ + a+-(S+S-)S+ , (6) 
S z Sj~S t = cti x .i z (Sj~ S[ )S i , (7) 
S+S*Sf=a?(S+Sf)S*. (8) 



Here, a\ x and a\ z are attached to NN correlations in 
the cry-plane and along the z direction, respectively, and 
a 2 is associated with longer ranged correlation functions. 
We obtain 

«-) 2 = [1 + 2a+-(C+ - + 2C+r, )](l - A 7q ) 

+A[1 + 4a+-((7 2 - 0i0 + 2Cf* li0 )](A - 7q ) 
+2a+-C+ - [A(4 7 2-l)-3 7q ] 
+4a+-C^ 0i0 [4 7 2-l-3A 7q ] 
R, 



(l + 2a+-C + - 2 )(l- Acos 9z ) 

+A(1 + 4a+-C z ^ 2 )(A - cosq z ) 
+2af~C^Q 1 (Acos2q z - cosq z ) 

+4a+-C^ 0A (cos2q z - Acos 9z ) 

2R Z 2a+-C+ fil (2- A 7q - Acosfe) 
+4a+-Cff 0il A(2A - 7q - cosq z ) 
+2a+~C+ ; i7q (Acos^ - 1) 
+4a+-C z f 0il7q (cos gz - A) 
+2atx C tofi cos^(A 7q - 1) 
+4a+-Cf 0i0 cos 9z ( 7q - A)] , 

«) 2 = 2(1 - 7q ) [l + 2af (C+ - + 2C+r i0 ) 
-2AagC+ - (l + 4 7q )" 
+R 2 z (l - cosq z ) [l + 2af C+ - 2 
-2Aa-C + - 1 (l + 2cos fc ) 

+8i? z [af C 1 + . 1 (2 - 7q - cosg z ) 
+Aaf|C^ l7q (cos^ - 1) 
+Aal z x C'l^ cos^( 7q - 1) . 



(9) 



(10) 



Note that in the special cases R z = and A = 1 the spec- 
tra reduce to the expressions given in Refs. [13] and [9], 
respectively. 
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The LRO in the correlators is reflected in our the- 
ory by the closure of the spectrum gap at Q = (w, n, it) as 
T approaches from above, so that liniT->T^ (Xq) = 
and Wq = at T < T^. Separating the condensation 

part Ce'* from [cf. Eq. (5)], the magnetization ra v 
is calculated as 



(11) 



Considering the uniform static longitudinal suscep- 
tibility Xo z = lim q ^ ^qV(< z ) 2 , the ratio of the 
anisotropic functions M^ z and (i^q Z ) 2 = c 2 xy {q x + q 2 ) + 
c z q z must be isotropic in the limit q — > 0. This yields 
the condition 

{Cz/Cxy) 2 — RzCofi,l/C 1,0,0 (12) 

with the squared spin-wave velocities 

c xy — 2 + a 2 Z {^2,o,o + ^1,1,0) ~ 5Aaf^C^ 

+2R z (a?C+ - 1 Aa^C+ - ) , (13) 

+4R z (arC+o,i - &(*llCo,o,\) ■ (14) 

Concerning the vertex parameters in our self- 
consistency scheme, three parameters are fixed by the 
sum rules' C* +~ = 1/2, C§* = 1/4, and by Eq. (12) 
for all T. To determine the free parameters taken as 
a tx~iz an( l a I Z j we need additional conditions. Let us 
consider the ground-state energy per site which we com- 
pose approximately, following Ref. [9], as e(A,R z ) = 
e(A, R z ) + e(A, 1) - e(A, 0). Here, e(A, R z ) denotes the 
ground-state energy of the 2D spatially anisotropic XXZ 
model [Eq. (1) without sum in y direction], where the 
values in the ID (R z = 0) and 2D cases (R z = 1) are 
taken from Ref. [22] and the exact data of Refs. [13,14], 
respectively. Since e(A,R z ) is known for R z = 0, 1 and 
A = 1 (taken from the Ising-expansion results by Affleck 
et al. [23]), we approximate e(A,i? z ) by the linear inter- 
polation e(A, R z ) = e(A, 0) + [e(A, 1) - e(A, 0)] [e(l, 1) - 
e(l,0)]- 1 [e(l,i? z )-e(l,0)]. At T = 0, we adjust af~ to 
e(A, R z ) and af to de/dA = 2C^ + R Z C^ Z 01 . 

To formulate conditions for and af z also at finite 
temperatures, we follow the reasonings of Refs. [21,10,12]. 
That means, we conjecture that the "vertex corrections" 
a ix(T) — 1 an( l o^'iT) — 1 have similar temperature de- 
pendences and vanish in the high-T limit. Correspond- 
ingly, as the simplest interpolation between high tem- 
peratures and T = we assume the ratio of two vertex 
corrections as temperature independent and fixed by the 
ground-state value, i.e., 



"(T)-l 



"SCO - 1 

(T)-l 



a 



(T)-l 



const. 



= const. 



(15) 
(16) 



To determine af~ (T) , as compared with af z z (T) fixed by 
the "isotropy condition" (12) resulting from Xo S we first 
note that an analogous condition cannot be derived from 
Xo~ = lim q ^ Mq ~) 2 j since both Mq 1 ' and o>+~ 

have non-zero q — > limits [cf. Eqs. (3) and (9)]. There- 
fore, for af~(T) we make the plausible ansatz assuming 
the ratio a" z (T)/a" x (T) as v independent, i.e., 



a 



lz 



CO af*(T) 



4x(t) 



(17) 



From the solution of the self-consistency equations the 
AFM correlation lengths above Tjy may be evaluated. 
They are obtained by the expansion of \q around Q, 
X^X^[l + {^ v ?{kl + kl) + {^fkl]^ withk = q-Q. 
The squared intraplane correlation lengths are given by 

(4V) 2 = -K")" 2 (is) 

X T + ^ a 2 (l^2,0,0 + ^2,0,0 + ^Ll.O + 2Cf Z 10 ) 

+a+" ((4A + f )C+ - + (8 + 3A)CT 0i0 ) 
+R Z [«+" (AC+ - + 2C- >0 ) + at~ (Cfci + 2Q,o,i) 
+A( a +-(C + - 1 + 2C ZZ 0A ) + a+-(C+^ + 2Cf* 0il ))]' 
-(AC+ - + 2C 1 - 0j0 )/M+-, 

x 'i + af (C+ - + 2C+r ;0 ) + HAaf*C+ - 

+2i? z (af z AC+ ; + 2aL z AC + - 1 + af C^) 
-2C+ - /M% . 
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